Ultra-high energy nuclei source in the direction to Virgo cluster. 
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The significant anisotropy in the arrival directions of the 69 events with energy E > 55 EeV 
detected by Pierre Auger collaboration [J is located in the 20-degree region centered near Cen A. 
Not only the 2-point, but also the 3-point and 4-point autocorrelation functions are completely 
saturated by this region. Besides there is an deficit of events in the direction of Virgo cluster. If one 
assumes that the excess around Cen A is due to heavy nuclei shifted from Virgo, one can explain 
20-degree scale of this anomaly. Also location of the highest energy event between the Cen A region 
and the Virgo cluster supports this idea. Magnitude and direction of the magnetic field is similar 
in this case to those expected for Galactic models. The existence of nuclei sources in the sky opens 
the road for a self-consistent description of Auger data. 



Introduction. Recently a very important step towards 
an understanding of Ulra-High Energy Cosmic Rays 
(UHECR) was done by Pierre Auger Collaboration [2J. 
For the first time, the large statistics of hybrid data al- 
lowed to study composition of UHECR not only using the 
average maximum depth of air-showers < X max >, but 
also the width of distribution, RM S(X max ). The main 
result of this study is that at energies around E > 10 EeV 
the composition changes from mixed to one compatible 
with heavy nuclei at E > 30 EeV Q . This result is con- 
sistent with observations of the average value of X max . 
Still it is a puzzle why the measurements of the HiRes ex- 
periment are consistent with a proton composition Q, a 
result that is in line with preliminary studies of the Tele- 
scope Array (TA) One can argue that either this is 
still a statistical fluctuation, or experimental studies of 
composition have unresolved systematic uncertainties or 
one starts to see a real North-South asymmetry in the 
UHECR composition. Measurements in the near future 
and detailed comparison of Auger to TA results will clar- 
ify this issue. 

There are several anisotropy signals present in the 
Auger data. A correlations between the arrival direc- 
tions of 27 highest energy events in the Pierre Auger Ob- 
servatory data and nearby AGN's from the 12th Veron 
catalog [5j was found in the prescribed test of Auger with 
the chance probability at 1 % level @. With increasing 
data set from 27 events in Q to 69 events in ref. [l| this 
correlation stayed at the 1 % level. Thus it significance 
strongly decreased. 

Also the autocorrelations in the Auger data up to Au- 
gust 2007 calculated with the two-point autocorrelation 
function showed an anisotropy with 1.6% probability by 
chance Q. This probability was reduced to 10% in the 
new Auger data with 69 events [l[. 

From other side, the excess of events in the 18° circle 
centered on the Cen A position is with 4% probability by 
chance as shown by Kolmogorov-Smirnov test [l|. This 
excess was already present in the data with 27 events, as 
was first mentioned in the ref. Q. 

Finally, there is a deficit of events towards the Virgo 



cluster. It is not very significant for a uniform distribu- 
tion of events on the sky, since Virgo is located at the edge 
of the Auger exposure. However, the situation is different 
for proton sources whose density follows the large scale 
structure. In this case, one expects a strong signal from 
this region, as was first discussed in the ref. [9j. Shift- 
ing the global energy scale by 20-30 % up as required by 
self-consistency of the muon data [Io[ will increase also 
the significance of the Virgo anomaly. 

The goal of this paper is to provide self-consistent de- 
scription of the Auger data within one model, both the 
anisotropy and the composition at the highest energies. 
Following the composition result we assume that most of 
69 UHECR observed at E > 55 EeV are heavy nuclei. 
In this case typical deflections of UHECR in the galactic 
magnetic field are about 50-100 degrees for E > 55 EeV 
nuclei [H|| . Thus one can not be surprised that there are 
no events coming from Virgo direction, since they are 
deflected away. 

Autocorrelations in the Auger data. First we study 
possible anisotropics in the Auger data in the standard 
way with autocorrelation functions. An autocorrelation 
signal at medium scales S ~ 25° was already presented 
in the analysis of combined data of experiments before 
Auger [13]. The two-point cumulative autocorrelation 
function as function of angle had a broad minimum for 
the first 27 Auger events with probability to be a fluctu- 
ation of background of the order of 10~ 2 ■ However for 
69 events this function shows an anisotropy with chance 
probability of 10% only fl. The existence of a cluster 
of events around Cen A mean that the anisotropy can 
remain visible in higher order autocorrelation functions, 
which arc more sensitive to clusters of events. 

To test this, we define a three-point cumulative auto- 
correlation function in analogy with the two-point func- 
tion. For given angle, we count all triplets of events which 
are separated by less then this angle. In practice, we cal- 
culate 3 angular distances between any three events and 
find the maximal one. In the same way one can define n- 
point correlation function with n points located at given 
angular distance from each other. Here due to computa- 
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tion reasons we stop on 4-point function. Then we sum 
up the number of n-point contributions in a given angle 
in the data and compare it to N = 10 5 Monte Carlo sim- 
ulated data sets with 69 events following Auger exposure. 
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FIG. 1: Probability that signal in the autocorrelation function 
is fluctuation of background as function of angle for first 27 
Auger events (thin lines) and complete up to date set of 69 
events with E > 55 EeV (thick lines). Blue dashed line is 
for two-point autocorrelation function, magenta dotted line is 
for three point autocorrelation function and red solid line for 
4-point autocorrelation function. 

Results for 69 Auger events up to December 2009 [l[ 
are compared to those with first 27 events, published in 
2007 in Fig. [TJ One can see that contrary to the two- 
point correlation function the three-point and four-point 
correlation functions are reduced only by the factor 2 
or 3 and still have a sharp minimum at 18-24 degree 
scale. Penalization in energy can not be calculated with- 
out knowledge of low-energy events with E < 55 EeV. 
However, one can estimate that penalty factor in anal- 
ogy with two-point correlation function, which is around 
100 in this case. This give 3% chance probability for 
3-point and 2% for 4-point autocorrelation functions. 

Since the angular scale in the 3 and 4-point autocor- 
relation functions is similar to the angular size of cluster 
of events near Cen A, we can check how significant is the 
contribution of this cluster to those functions. In order 
to do this we cut a 20-degree region around Cen A both 
from data and from MC simulations and calculate again 
the three-point and 4-point autocorrelation function. 

The result is presented in Fig. [2] It is seen from this 
figure that ALL signal in the 2, 3 and 4-point autocor- 
relation functions come from this single cluster of events 
near Cen A. Moreover, thin lines show that data without 
this cluster of events is too isotropic on the scale of clus- 
ter around 20 degrees as compared to random MC (up 
to few percent level). 

Thus the autocorrelation analysis of 69 Auger events 
shows that all existing anisotropy in data itself (if not by 
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FIG. 2: Probability that signal in the autocorrelation function 
is fluctuation of background as function of angle for set of 69 
events (thick lines) and set of 56 events excluding 20-degree 
region around Centaurus A (thin lines). Line code same as in 
Fig. ffl 



chance) is due to a single cluster of events near Cen A, 
which we call below Cen A anomaly. 

Nuclei source in direction to Virgo cluster and Cen A 
anomaly. Now let us recall that there is another anomaly 
in the Auger arrival directions. There are no events in the 
20-degree region around the Virgo cluster. This can be 
a problem if proton sources are following the matter dis- 
tribution at the Large Scale Structure [9j. However, the 
expected number of events from this region depends on 
many assumptions about sources and their distribution. 
Under the assumptions of ref. [§] one expect 15 events 
out of 69 coming from this region (which corresponds to 
original 6 out of 27). In the other extreme of models in 
ref. [l| one expect just few events. Detailed discussion 
of this subject is beyond scope of present study, since in 
the case of nuclear sources this anomaly would not be a 
problem. Nuclei coming from the Virgo region would be 
shifted to some other part of the sky. Below we assume 
that the events around Cen A are those which come from 
the Virgo cluster. Then due to the difference in exposure 
in 2.5 times between these two regions one expect that 
13 nuclei in Cen A would correspond to 5 proton events 
in direction to Virgo, which is in the range of models 
discussed above. 

Such an assumption can be tested qualitatively in sev- 
eral ways. First, one has to explain the 60-degree shift 
from Virgo to Cen A and the 20-degree spread of events in 
the cluster. Second, deflections can have partial informa- 
tion on 5 ~ Q/E for very different energies, despite this 
dependence can be non-trivial [111 ]. Third, the required 
parameters of Galactic magnetic field should obey exist- 
ing constraints. Let us show that our assumption pass 
all above tests. 

Deflections of UHECR nuclei with E > 60 EeV in the 
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FIG. 3: Sky map with 69 Auger events in the in the galactic 
coordinates. Stars show events outside of 20-degree region 
around Centaurus A, squares inside. Red color for events 
with E > 100 EeV, magenta for 70 EeV < E < 100 EeV and 
blue for 55 EeV < E < 70 EeV. Positions of Cen A and Virgo 
are shown with green points. Blue arrow show direction of 
deflection from center of Virgo towards highest energy event 
due to Galactic Magnetic Field. 



regular galactic magnetic field were studied recently in 
Ref. " 
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Typical deflections arc expected about 50-100 
degrees depending on the source position in the sky. Thus 
60-degree deflections arc natural in this case. Here let us 
note a possible explanations of 20 degree spread of events 
around Cen A. It can be the sum of several effects, in- 
cluding deflection in the turbulent magnetic field, the dis- 
tribution of sources in the Virgo region and contribution 
of large deflections of nuclei in the global magnetic field 
of Virgo cluster In the last case even if galaxy M87 
is the single source of nuclei in Virgo region, its signal 
would be spread over 10-degree scale - size of Virgo clus- 
ter on sky, before it reach magnetic field of our Galaxy. 
Besides extragalactic field between Virgo and our Galaxy 
can distort the signal as well. 

Second, due to the fact that the deflection in regular 
field (60 degrees) is bigger than the spread of events (20 
degrees) one expect that the events around Cen A are 
ordered according to their energy. Events with higher 
energy should be deflected less. In order to test this we 
divide all our events with E > 55 EeV in three bins, 
E > 100 EeV, 70 EcV< E < 100 EeV and 55 EcV< E < 
70 EeV. Results are presented in the Fig. [3l There are 4 
events in 70-100 EeV bin and remaining 9 events in lowest 
energy bin. The distribution of events within cluster is 
unclear. However, the highest energy event with E = 140 
EeV is located between Cen A and Virgo in agreement 
with our assumption that it originate from Virgo. Let us 
note that nuclei with charge Z should not be deflected 
strictly as 5 ~ Z/E, which is true only in case of protons 
(see examples in ref. [HI). 



Finally, let us discuss deflections in the Galactic and 
Extragalactic Magnetic fields. If the picture above is cor- 
rect, one can estimate the value of the magnetic field 
responsible for the deflections. For this we will use the 
highest energy event. Its energy is E = 1.4 x 10 20 eV 
and if it was emitted in the center of Virgo cluster the 
deflection was S = 34°. 

For deflection by a regular field we have the relation 



34° = Segmf + 1° • Z- 



B 



1.4 • 10 20 eV 



kpc 3-10" 6 G 



E 



, (1) 



where Segmf is part of deflection due to extragalactic 
magnetic field, d is the distance which particle travelled 
in a Galactic Magnetic Field (GMF) of strength B. E 
and Z are the particle energy and charge. 

In general, Eq.((T]) can have two solutions, depending 
on value of Segmf- If Segmf ~ 30° for protons as in 
ref. [l4[, one can describe data with proton deflections 
Q = 1 dominated by extragalactic field. However, if 
Segmf 30° as in ref. [TBI one can have a solution with 
heavy nuclei Q < 26. Auger data on composition, show 
that the composition becomes heavy at high energies Q 
however one can not exclude protons from consideration 
before confirmation of this result. 

Let us note that not only the value but also the re- 
quired direction of GMF is good, since deflections are 
almost perpendicular to the galactic plane. As it is well 
known, z-componcnt of GMF is B z = 0.2/iG, which is 
much smaller than B x ^ y ~ 2/iG. This leads to UHECR 
deflections almost perpendicular to the galactic plane. 
Most of existing GMF models obey this condition. And 
Auger events in Fig. [3] fulfill such condition. 

Here let us add few words about UHECR sources. One 
can expect that only exceptional objects can accelerate 
particles to highest energies [l6|. In this respect Cen A 
galaxy is not a very powerful object 17j. To the con- 



trary, M87 galaxy, as center of the local Large Scale 
Structure and Virgo cluster was suggested as UHECR 
source long ago [18j. Probably we start to see its sig- 
nal in UHECR nuclei. However, realistic discrimination 
of source(s) in the Virgo region will require much larger 
statistics as compared to present one. Note that UHECR 
flux of Auger events in Cen A region is of the same order 
as the gamma-ray flux from M 87 measured by HEGRA 
F ~ 10" 12 erg/ cm 2 /s 0. 

Another important constraint on the properties of the 
nuclei source responsible for the Cen A excess was found 
in Ref. [20[ ■ Protons with energy E = 2 — 3 EeV have the 
same rigidity as E = 60 EeV iron nuclei. Such protons 
would be therefore deflected in same way as heavy nu- 
clei, independent of the details of the unknown magnetic 
fields between the source and the Earth. At E = 2 — 3 
EeV, a la fluctuation of the Auger background corre- 
sponds to about 100 events [2(J. Thus if the heavy nuclei 
source also accelerates protons, they can be visible above 
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the background. Since the number of events at the high- 
est energies is around 10, the constraint of Ref. [2(J does 
not allow a soft spectrum of accelerated particles and a 
very large proton to nuclei ratio, which combined should 
not lead to a proton flux (at E = 2 — 3 EeV) exceeding 
the UHECR nuclei flux (at E = 60 EeV) by a factor 10. 
Since al/E 2 spectrum already gives a factor 10, this con- 
dition requires either a harder spectrum or not more than 
one to one ratio of accelerated protons to nuclei at the 
source. Note that the magnetic field of the cluster can 
significantly change the original power law spectrum in a 
given direction [l3j. Another option can be acceleration 
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with a very hard spectrum, as in the model of Ref. 
Thus the results of Ref. [20| show that the confirmation 
of a nuclei source in Virgo would inform us also about the 
acceleration mechanism and would constrain the compo- 
sition of accelerated particles in this source. 

Conclusions. We found that the autocorrelation signal 
in recent Auger data [l[ still show some anisotropy in the 
3-point and 4-point autocorrelation functions on a few 
percent level. We found that this signal is completely 
due to 20-degree scale cluster of events located near Cen 
A and there is no anisotropy remaining after removing 
this cluster of events from data. 

Another anomaly in the UHECR arrival directions in 
Auger data located near Virgo. There are no events 
with E > 55 EeV coming within 20-degree region of this 
galaxy cluster, which is in conflict with expectations if 
proton sources are following Large Scale Structure @ . 

Since Auger see heavy nuclei dominated composition at 
highest energies Q, we suggest that both anomalies can 
be explained by a heavy nuclei source located in Virgo 
region, with nuclei shifted from it towards Cen A by the 
Galactic Magnetic Field. 

In this case the highest energy event located between 
Virgo and Cen A belongs to this source. Both the 60- 
degree shift and the 20-degree spread of events can be 
explained by deflections in the galactic magnetic field. 
Required fiG magnitude and direction of GMF with 
\BJ < \B. 



are comparable to expectations for GMF 
models. 

In the future Auger South can test the statistical signif- 
icance of the cluster near Cen A by a blind test. Another 
key point would be the collection of additional data on 
composition. In case a nuclei-dominated composition will 
be confirmed, one can test the hypothesis of nuclei from 
Virgo, presented here by a collection of a large fraction 
of E > 10 20 eV events, some of which have to be located 
between the Cen A and Virgo regions. The required large 
statistics will be possible for future JEM-EUSO space ex- 
periment. Also, if similar anomalies will be found by the 
Telescope Array in the North Hemisphere, Auger North 
will be able to study them. 
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